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' Abstract. We study the multi-wavelength properties of a set of 171 Lya emitting candidates at redshift z — 2.25 found 

in the COSMOS field, with the aim of understanding the underlying stellar populations in the galaxies. We especially seek 
to understand what the dust contents, ages and stellar masses of the galaxies are, and how they relate to similar properties 
of Lya emitters at other redshifts. The candidates here are shown to have different properties from those of Lya emitters 
found at higher redshift, by fitting the spectral energy distributions (SEDs) using a Monte-Carlo Markov-Chain technique and 
including nebular emission in the spectra. The stellar masses, and possibly the dust contents, are higher, with stellar masses 
| in the range log M* = 8.5 — 11.0 M© and Av = 0.0 — 2.5 mag. Young population ages are well constrained, but the ages 

f"^ . of older populations are typically unconstrained. In 15% of the galaxies only a single, young population of stars is observed. 

We show that the Lya fluxes of the best fit galaxies are correlated with their dust properties, with higher dust extinction in 
Lya faint galaxies. Testing for whether results derived from a light-weighted stack of objects correlate to those found when 
■ fitting individual objects we see that stellar masses are robust to stacking, but ages and especially dust extinctions are derived 

incorrectly from stacks. We conclude that the stellar properties of Lya emitters at z — 2.25 are different from those at higher 
redshift and that they are diverse. Lya selection appears to be tracing systematically different galaxies at different redshifts. 
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1 . Introduction group of high redshift galaxies are found based on their red 

colours (DRGs, EROs; Franx et al. 2003; Daddi et al. 2004). 

The high redshift Universe hosts a zoo of various types of Unlike ^ LBGf , and ag expected from their colourSi ^ 

galaxies, selected with many different methods. Based on num- , ■ e . ■ , , . c , u 

fe 3 galaxies are often massive and dusty (Forster-Schreiber et 

bers of detections, the by far most abundant group of star- ^ 2QQ4 . Papowich et al 2006) Qther types of high redshift 

forming galaxies with redshifts z > 3 are the Lyman-Break , axies star _ bursti laxies such as those detected 

Galaxies (LBGs; Steidel et al. 1996, 1999, Madau et al. 1996, wift sub . mm arrayfj (SMGs; Mson et ^ ^ ^ Blain et 

Papovich et al. 2001, Giavalisco 2002, Ouchi et al. 2004), ^ m% chapman et al 2005) and smaller galaxies detected 

found by searching for the Lyman-break in the spectra of galax- ■ u . . /D1 . , lriri o 

3 & 3 when gamma-ray bursts occur in them (Bloom et al. 1998, 

ies using broad-band imaging. These galaxies often turn out to Fynbo ^ ^ 2QQ ^ christensen et al 2004) Qne selection 

be moderately star-forming, medium mass galaxies (Shapley et method ^ has become increasingly popular is that of narrow . 

al. 2005; Reddy et al. 2006; Verma et al. 2007). Another large band imaging of ^ highly redshifted Lya line in star . forming 
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* Based on observations carried out at the European Southern 

/~>u , , CCA > i TT ^ XT „„. , m ,„, D , between redshifts z = 2 — 5 (Gronwall et al. 2007, Venemans 

Observatory (ESO) under prog. ID No. 084.A-0318(B). v 



galaxies (e.g. M0ller & Warren 1998, Cowie & Hu 1998). The 
largest samples to date include several hundred Lya emitters 
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et al. 2007, Nilsson et al. 2007, Finkelstein et al. 2007, Ouchi 
et al. 2008, Grove et al. 2009, Guaita et al. 2010). Slowly, stud- 
ies of the actual properties of the type of galaxies selected with 
this method, most notably through SED fitting, are getting un- 
der way. 

There is a large spread in the results of previously published 
work on the topic of Lya emitter properties and much work is 
still needed to determine what type of galaxies are selected with 
the Lya technique, and where they belong in the high redshift 
galaxy zoo. To date, eleven publications have presented the 
SED fitted properties of high redshift Lya emitters (Gawiser et 
al. 2006, 2007; Finkelstein et al. 2007, 2009a; Lai et al. 2007, 
2008; Nilsson et al. 2007; Pirzkal et al. 2007; Yuma et al. 2010; 
Ono et al. 2010a, 2010b), and two in the local Universe (Hayes 
et al. 2007; Finkelstein et al. 2009b). Lya emitters are selected 
based on their properties in the UV and have often been found 
to be extremely faint at longer wavelengths. Since the UV light 
is dominated by young stars, also contributing to the Lya emis- 
sion, little is thus known about the properties of any underly- 
ing populations of stars. For galaxies above redshift two, it is 
thus crucial to have detections of the galaxies in the near- and 
mid-IR, corresponding to restframe optical and near-IR wave- 
lengths. Few Lya emitters have been detected in these wave- 
length regimes, with only a handful to date at z ~ 3 (Lai 
et al. 2008), and at z = 4-6 (Lai et al. 2007, Pirzkal et 
al. 2007, Finkelstein et al. 2008, 2009a, Yuma et al. 2010). 
These studies have shown IR-detected Lya emitters to be in 
general more massive, and in some cases more dusty than their 
non-IR-detected counter-parts. It is not yet clear whether there 
are indeed two (or multiple) populations of Lya emitters (e.g. 
faint-small-dust-free vs. bright-massive-dusty) or if there is a 
smooth transition of objects with various masses and proper- 
ties. 

The surveys to date have concentrated on redshifts z > 3. 
In an effort to get a better handle on the SED of the selected 
galaxies, we here present the SED study of a sample of Lya 
emitters at lower redshift, z = 2.25, first presented in Nilsson 
et al. (2009a). The advantages of the lower redshift include the 
possibility of observing a larger range of the restframe spectra 
with the same instruments, and increased signal-to-noise in the 
detections due to the smaller luminosity distance. Indications 
of evolution in the properties of this type of galaxy from red- 
shift z ~ 3 to z ~ 2 were reported in Nilsson et al. (2009a); in- 
cluding typically redder colours, smaller equivalent widths and 
a larger AGN content in the sample. In this paper we investigate 
the stellar properties of the LAEs further by fitting their SEDs 
with theoretical models. The question we seek to answer is how 
these properties have evolved from higher to lower redshift. 

This paper is organised as follows: section [2] describes the 
sample of Lya emitters that are being studied, and includes an 
analysis of the Spitzer detections of the galaxies. In section [3] 
we outline the method used for the SED fitting, and the results 
of the fitting are reported in section [4] A discussion of the red- 
shift evolution of LAE stellar properties is found in section [5] 
We end with a general discussion in section|6]and a conclusion 
in section Q 



Throughout this paper, we assume a cosmology with Hq = 
72 km s _1 Mpc~\ Q m = 0.3 and Q\ = 0.7. Magnitudes are 
given in the AB system. 



2. Data 

2.1. Sample 

The sample studied in this paper is nearly identical to the sam- 
ple in Nilsson et al. (2009a). In brief, this sample was found 
through narrow-band imaging with the ES02.2m/MPG tele- 
scope on La Silla, Chile, using the Wide-Field Imager (WFI; 
Baade et al. 1999). A central section of the COSMOS field 
(RA. 10 h 00 m 27 s , Dec 02°12'22'/7, J2000) was observed with 
the filter N396/12, observing Lya with z = 2.206 - 2.312. The 
data reduction and method for selecting emission line objects 
is explained in detail in Nilsson et al. (2009a). In the survey, 
187 emission-line objects were found. Of these, 17 were de- 
tected in the public GALEX data available and were flagged 
as potential [Oil] emitters, although this sub-sample will also 
contain Lya emitters with high escape fraction in the restframe 
UV. The photometry of these candidates in the optical and near- 
infrared data available in COSMOS is presented in Nilsson et 
al. (2009a). 

In January/February 2010, multi-object spectroscopy with 
the VIMOS instrument on VLT UT3 (programme ID 084.A- 
0318(B)) was carried out on 152 candidates. The results of 
this spectroscopic campaign will be presented in a forthcom- 
ing publication, but we will in the analysis here exclude those 
candidates which are clearly confirmed to not be Lya emit- 
ters in the spectroscopy. Conversely, a number of less likely 
(named "unlikely" in Nilsson et al. 2009a) candidates were 
put on slits. Those that showed a Lya detection in the spec- 
troscopy will be added to the analysis here (a total of six). 
The results of the Nilsson et al. (2009a) paper do not change 
significantly from the exclusion/inclusion of the spectroscopic 
rejections/confirmations. In Table Q] we present the different 
sub-samples of the candidates in the various steps. In the fi- 
nal sample that will be used throughout the paper here-after, 
only candidates that are either not rejected in the spectroscopy, 
or did not get observed spectroscopically, are included. With 
"not rejected" we take to mean those candidates that were spec- 
troscopically observed and that were left after rejecting galax- 
ies with multiple spectral lines inconsistent with the Lya red- 
shift, or which did not have significant line emission on a de- 
tected continuum. We further add a few more AGN by study- 
ing the Chandra data in the COSMOS field (see sec. 12.3b . re- 
ducing the number of AGN-free LAE objects. The final sam- 
ple thus consists of 171 LAEs, of which 141 are AGN- and 
GALEX-undetected. Note that at least 8 GALEX-detected non- 
AGN LAEs are confirmed as being true Lya-emitters. These 
represent a unique sub-sample of high-transmission objects at 
high redshift (in agreement with M0ller & Jakobsen 1990), and 
would be worthy a special follow-up. 
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Table 1. Definitions of all sub-sets of LAE candidates. 
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Sample 


Total 


LAE best 


GALEX only 


AGN only 


GALEX + AGN 


Original 


187 


156 


12 


14 


5 


Original (C-list) 


61 










Spec. obs. 


118 


96 


11 


8 


3 


Spec. obs. (C-list) 


34 










Included f. spec. 


96 


83 


8 


4 


1 


Included f. spec. (C-list) 


6 


6 











Total after spec. 


171 


149 


9 


10 


3 


Final (after new AGN) 


171 


141 


7 


18 


5 


Fitted 


58 


54 


4 






Good Xr 


40 


38 


2 







Each line defines a further sub-sample of the original Nilsson et al. (2009a) sample (first lines; divided into the originally published sample 
and the C-list, or "unlikely" candidates from that paper), divided into the total sample, the sample without GALEX and AGN detections, 
and the GALEX only, AGN only and combined GALEX and AGN detected LAEs. The "unlikely" candidates were not analysed for AGN 
and GALEX detections prior to spectroscopic confirmation. The second to fifth rows show the number of candidates put on slits, and those 
that were not excluded. The row titled "Total after spec." indicates the resulting sample after the spectroscopy, summing those objects 
that are accepted by the spectroscopy with those without spectroscopic data, and the following row indicates the result of the new AGN 
detections (sec. 12.3b . The final rows show how many of each sub-sample were fitted for and which had good Xr m the fit (sec. l4.3t . 



2.2. Updated photometry 

Since the publication of the original sample, several new im- 
ages have become available from the COSMOS consortium, 
including J images from UKIRT, K s images from CFHT and 
X-ray data from the Chandra observatory (see sec. [23). The 
depths of the infrared images are 23.7 magnitudes (5er) in both 
bands (Ilbert et al. 2009). We prepared and extracted fluxes of 
all the objects in the catalogue in a similar way as described in 
Nilsson et al. (2009a). The number of detections was 43 and 
86 of the 171 candidates respectively in the J and K s bands. 
One of the spectroscopically confirmed, less likely candidates 
was detected in the K s band. These results were incorporated 
into our main photometric catalogue. In the case that an object 
was already detected in the previously available, but shallower, 
KPNO K s band, this magnitude was replaced with the CFHT 
result for consistency. 

For all optical and near-IR bands used in the fitting (Bj, 
Vj, i, z, J, and Ks), the curves of growth of 25 stars in the 
field were measured and aperture corrections from 3" diame- 
ter apertures to "infinity" calculated. These were applied to all 
magnitudes before the SED fitting. We also subtracted the ob- 
served Lya flux (as calculated from the narrow-band observa- 
tions, Mag3ggo) from the Bj band, as this line is exceptionally 
difficult to fit. 

2.3. Chandra detections 

Elvis et al. (2009) presented the available Chandra X-ray im- 
ages, with flux limits of 1.9, 7.3 and 5.7xl0~ 16 erg s _1 cm~ 2 
in the soft (0.5 - 2 keV), hard (2 - 10 keV) and total (0.5 - 
10 keV) bands respectively. To search for detections, aperture 
photometry with 5" diameter was made on the positions of the 
candidates in the 0.5 — 10 keV band. For detections at the So- 
level, the photometry was repeated in the soft and hard bands. 



This revealed 19 detections, of which five are only detected in 
the total band, five/two are detected in the total and soft/hard 
bands. Seven candidates are detected in all bands. Based on the 
1761 sources found in the whole 0.5 deg 2 COSMOS survey 
by Elvis et al. (2009), the random chance of an object in one 
aperture is 0.0053, or 0.9 in the whole sample of 171 objects. 
All LAEs detected with Chandra (or previously with XMM or 
VLA, see Nilsson et al. 2009a) are considered AGN due to their 
bright X-ray (or radio) fluxes. The faintest limit in the X-ray 
data is that of the soft band from Chandra, and this limit corre- 
sponds to a SFR of ~ 1500 M© yr _1 according to the conver- 
sion factor of Ranalli et al. (2003). Similarly, the faintest radio 
detection corresponds to a SFR of ~ 950 M Q yr" 1 according 
to the conversion of Bell (2003). These limits are much higher 
than those derived from the Lya line in these galaxies, with the 
conclusion that all X-ray /radio detections must be AGN-driven. 
Their R band-to-X-ray fluxes are in the range of what is com- 
monly found for AGN (Schmidt et al. 1998, Ffornschemeier et 
al. 2001). Hence we now have a total of 23 LAE AGN, 10 of 
which are only detected with Chandra, three are only detected 
in radio emission with VLA, one is only detected with XMM, 
and 9 that are detected by both Chandra and VLA/XMM. The 
new AGN detections are LAE_COSMOS# 14, 20, 52, 67, 76, 
120, 130, 134, 151 and 175. 

Thus the confirmed AGN contamination in the non- 
GALEX-detected sample increases from 9% (before Chandra 
detections but after spectroscopic rejections) to at least 11%, 
i.e. 18/159 (see also Table [TJ, a much larger percentage than 
in high redshift samples (Wang et al. 2004, Gawiser et al. 
2007). Lower redshift samples of Lya emitters at z = 0.3 and 
z = 2.2 have also been shown to contain large AGN frac- 
tions, of e.g. ~ 43% at z = 0.3 (Finkelstein et al. 2009c) 
and ~ 75% at z = 2.2 (Bongiovanni et al. 2010) indicating 
a clear redshift evolution, although the result of Bongiovanni 
et al. (2010) is based on detecting AGN with a method dif- 
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Table 2. S-COSMOS sensitivities from Sanders et al. (2007). 



Band 


CWL 


FWHM 


Sensitivity 


Chi 


3.58 


0.75 


0.0009 


Ch2 


4.50 


1.02 


0.0017 


Ch3 


5.80 


1.43 


0.0113 


Ch4 


8.00 


2.91 


0.0146 


MIPS (deep) 


24.0 


4.7 


0.071 


MIPS (shallow) 


24.0 


4.7 


0.15 



Central wavelengths and full width half maximums are given in 
/j,m. Sensitivities are 5cr, and given in mjy. Only MIPS 24 /im 
values are given as our candidates have no detections in the longer 
wavelength bands. 



ferent from other groups, which may influence the reliabil- 
ity of this fraction in comparison with other results. The in- 
crease in AGN fraction was discussed already in Nilsson et 
al. (2009a), and we here repeat their test of studying only a 
sub-sample of the LAEs that are bright enough to have X-ray 
detections if they are AGN. Using the same quasar compos- 
ite spectrum as in Nilsson et al. (2009a), the limiting magni- 
tude in Lya is 24.79 based on the new X-ray detection limit of 
the Chandra observations. Of the final sample in this publica- 
tion, 53 LAEs fulfill this criterion. Based on these LAEs only, 
the AGN fractions are 22 ± 7% and 26 ± 8% respectively, if 
the GALEX detected objects are excluded or included. Even 
though these fractions are lower than in Nilsson et al. (2009a), 
they are still ~ 2.5 times that what was reported at z = 3 in 
Ouchi et al. (2008). For a longer discussion on AGN fractions, 
see Nilsson & M0ller (201 1). The LAE AGN have total X-ray 
fluxes between log Lo.5-iokcV = 43.3 — 45.5 erg s _1 , if at 
z = 2.25. 



2.4. Spitzer detections 

The COSMOS field is covered by Spitzer in the S-COSMOS 
survey (Sanders et al. 2007). The public data includes the four 
IRAC channels (3.6, 4.5, 5.8 and 8.0 ^m), and the three MIPS 
bands (24, 70 and 160 /im). We will not discuss the 70 and 
160 /im bands further as no detections were made in those 
bands. In the MIPS 24 /.tin band, both a full and a deep survey 
have been published, with the deep survey covering roughly 
33% of the original narrow-band image used for selection of 
candidates. The sensitivities of the Spitzer data can be found in 
Table |2] We searched the public catalog (Sanders et al. 2007, 
v. "G02", released May 2007) for sources associated with our 
candidates. In the IRAC bands, counterparts were searched for 
within a circle with 5 IRAC pixels (3") radius around each can- 
didate, and in the MIPS bands within a circle with 4 MIPS pix- 
els (4.8") radius. In total [56, 41, 25, 15, 20] candidate counter- 
parts were found in the [Chi, Ch2, Ch3, Ch4, MIPS-24] bands, 
see also Table [3] that includes all sub-sets of the candidates. 
The random probability of finding an object within our search 
radii are 0.0096 in Chi, resulting in 1.6 random false detec- 
tions in the whole sample. After visually inspecting the associ- 



ations, two were excluded as random false detections and the 
list of Spitzer detected sources should, after removal of these 
two sources (detected in Chi and Ch2), be robust. The detec- 
tion rate of our candidates is 33% in Chi and 12% in MIPS-24, 
which is a higher rate than found in surveys for LAEs at red- 
shifts three and beyond (Lai et al. 2007, 2008) but in agreement 
with surveys at z ~ 2.3 (Colbert et al. 2006). The MIPS-24 de- 
tected candidates are by definition also ultra-luminous infrared 
galaxies (ULIRGs) due to their bright infrared fluxes. For a fur- 
ther discussion on these LAE ULIRGs and their evolution with 
redshift, see Nilsson & M0ller (2009). 

To test the public catalogues we performed aperture pho- 
tometry on a number of objects in the catalogue. The catalogue 
fluxes were extracted in l'/4 radius apertures (to be compared 
to l'/5 radius apertures for the optical broad-band images), and 
aperture corrected to give the total flux of the object accord- 
ing to the aperture corrections listed in the release notes of the 
catalogue. The fluxes measured in the image were in reason- 
able agreement with the catalogue values, however, the errors 
measured were in the cases of Chi - Ch3 roughly two times 
those quoted in the catalogue. For Ch4 and MIPS the errors 
were also in agreement. Thus we multiply the published errors 
on our candidates in the Chi - Ch3 bands by two for the fol- 
lowing analysis. There was a slight tendency of over-estimation 
of the fluxes in Chi, however, it was not significant enough to 
warrant a correction of the catalogue fluxes. It should be noted 
though, that the Chi fluxes may be over-estimated by up to 5%. 
In the subsequent SED fitting, an extra error of 10% of the flux 
of each object was included, in order to account for systematic 
errors in the data reduction (see also Muzzin et al. 2009). 

3. SED fitting method 

The method used to fit the SEDs of the LAEs is an updated ver- 
sion of the code used in Nilsson et al. (2007), called NisseFit. 
It is based on the stellar populations catalogue GALAXEV 
(Bruzual & Chariot 2003) and the fitting is done using a Monte- 
Carlo Markov-Chain (MCMC) algorithm. The algorithm ex- 
plores a multi-dimensional parameter space by stepping in a 
random fashion from one point Pi in this space to another Pi+i . 
The step is random in the sense that the projected distance from 
point Pi to Pi+i along one of the parameter axes in the param- 
eter space is a randomly chosen value ranging from to the 
maximum step size, AP max . The direction of the step along an 
axis is also random, with a probability of 50% that the direc- 
tion is along the positive or negative direction respectively. We 
impose boundary conditions for all parameters and ensure that 
no step leads to a value outside these boundaries. The param- 
eters and the boundary values we impose on them are listed 
in Table |4] Once a new point Pi+i is determined, the current 
model is advanced to the new point with a probability Pi^i+i 
which is given by the ratio of the two model probabilities given 
the data D, 

p{Pi+i\D) 



Pi^i+i 



p(P\D) 



Max(exp [( X 2 i+1 - x?)/T] , 1) 



(1) 

i+1 and xl are tne reduced \ 2 values for the two 
points in the parameter space and T is a "temperature" param- 



where xl 
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Table 3. Spitzer detections for all sub-sets of LAE candidates. 



Band 


Total sample 


LAE best 


GALEX only 


AGN only 


GALEX + AGN 


Chi 


56 (33%) 


36 (26%) 


2 (29%) 


13 (68%) 


5 (100%) 


Ch2 


41 (24%) 


24 (17%) 


(0%) 


12 (63%) 


5 (100%) 


Ch3 


25 (15%) 


8 (6%) 


(0%) 


12 (63%) 


5 (100%) 


Ch4 


15 (9%) 


5 (4%) 


(0%) 


7 (37%) 


3 (60%) 


MIPS 


20 (12%) 


8 (6%) 


(0%) 


7 (37%) 


5 (100%) 



For sample definitions, see Table[TJ Percentages of each sample are given in parentheses. 



eter which can be used to adjust the convergence. In our case 
T = 1. If the model is not advanced to the new point, the pa- 
rameters remain Pi, so Pi+i = P{. This algorithm has several 
advantages over more common maximization algorithms (such 
as steepest descend, amoeba, etc.) in that it allows a full ex- 
ploration of the parameter space. The result of our algorithm 
after N iterations is a list of points Pi,i = 0..N that, if N is 
sufficiently large, represents a good sampling of the underlying 
probability distribution function of the parameters. 

Please note carefully that, even though we refer to the quan- 
tity \r defined in Eq.[3]as the "reduced chi-squared" here and 
throughout this paper, the underlying modelling errors are not 
known to sufficient accuracy, compared to the relatively high 
measurement accuracy of the data, to allow the usual interpre- 
tation of a reduced chi-square value. Given the large number of 
photometric bands and the small measurement errors in many 
of these, no stellar population models are expected to give good 
fits to all the data in an absolute sense, and hence the \f- val- 
ues are not actually expected to be close to one. The term has 
been chosen here merely to reflect the fact that the "goodness 
of fit" is calculated in the same way as a reduced chi-squared. 
However, it should be interpreted as a relative merit value to 
compare different fits rather than an absolutely measurement 
of fit quality. 

3. 1. Improvement on models 

Recent publications show that nebular emission lines can play 
an extremely important role in fitting SEDs (Zackrisson et 
al. 2008, Schaerer & de Barros 2009, Raiter et al. 2010). In 
order to properly take into account this potentially large effect 
of gas emission in young, star-bursting galaxies, an add-on to 
the GALAXEV code was created. The nebular gas contributes 
with a continuum component and an emission-line component. 
The continuum addition was calculated from the Starburst99 
models (Leitherer et al. 1999) with a Salpeter IMF with slope 
a = 2.35 in the mass range 1 - 100 M© (the Salpeter IMF 
is also used in the GALAXEV models). The emission-line 
strengths were calculated using the MAPPINGS photoionisa- 
tion code (Kewley et al. 2011, in prep.). Both nebular contin- 
uum and emission-line strengths were calculated on a grid of 
ages (for the continuum in 36 steps between 1 and 100 Myrs 
and for the emission lines in six steps between 1 and 20 Myrs) 
and metallicities (same as available in GALAXEV, see below) 
and for each fit, the nearest values for these parameters were 




10 4 10 5 
Obs. wavelength (A) 



Fig. 1. Illustration of the added nebular emission. GALAXEV 
spectra at z = 2.25 with ages 1, 5, 20, and 100 Myrs are 
shown with nebular emission lines (red) and without (black). 
The spectra are scaled to arbitrary units to illustrate the effect 
of age on the nebular emission. The filtercurves of the filters 
used in the fitting (Bj, Vj, i + , z + , J, K s and Spitzer Chl-4) are 
shown at the bottom of the figure. The nebular emission lines 
are added as delta functions at the respective wavelength of the 
emission lines, and the step-like appearance of the 1 Myr spec- 
trum is due to the nebular continuum emission. At ages larger 
than 20 Myrs, the effect of the nebular add-on becomes negli- 
gible. 



used. The emission-lines were scaled with the Ha flux and the 
number of ionising photons in the galaxy spectra produced by 
GALAXEV. For an illustration of the added nebular emission, 
see Fig.Q] In practice, the emission line strengths as well as the 
continuum contribution become negligible at ages > 20 Myrs. 

In this publication, we choose to fit only two single stel- 
lar populations (SSP) models, as the nebular emission add- 
on becomes too complicated with any other star formation 
history. For all models, dust is added according to Calzetti 
et al. (2000). The parameter space allowed is detailed in 
Table |4] For the metallicity, the steps allowed are Z/Zq = 
0.005,0.02,0.2,0.4, 1.0. Only fluxes in the bands from Bj to 
the 8.0/im bands are used, as there are large uncertainties in 
the models at wavelengths below the Lya line, and at rest- 
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Table 4. Parameter space allowed by SED fitting code. 



100 





Min 


Max 


Remark 






Old pop. age 


OGyr 


2.7 Gyr 




CL 

o 


80 - 


Young pop. age 


OGyr 


0.1 Gyr 


< Old age 


Q_ 




Mass frac. in young pop. 
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frame mid-IR wavelengths. The observed Lya flux, as calcu- 
lated from the narrow-band photometry, was subtracted from 
the Bj band before fitting to avoid uncertainties in the fitting 
of the Lya line itself. Correspondingly, the Lya line was also 
subtracted from the nebular emission add-on. 

When the spectrum of a very young SSP is added to an 
older SSP spectrum, the older spectrum will need to comprise 
a very large mass fraction in order to be seen, as the young SSP 
will dominate the total spectrum. To test at what mass frac- 
tions the old population is actually seen, several young and old 
SSPs were generated and co-added with different mass frac- 
tions. To find an observability criterium for the older popu- 
lation, the colour in Vj— IRAC Chi was calculated. This was 
compared to the typical la error bar in the photometry in these 
two bands. We calculated the colour significance according to: 



col,i 9 nif.(mf) = ^ = Chl ^ - Z ^W=i-° 

(2) 

where mf is the mass fraction. The number in the denominator 
is the typical error on the colour based on the data at hand. If 
this colour significance is more than la, the old population is 
considered to be observed. A plot of these limits can be found 
in Fig. |2] Based on this figure, all objects with more than 90% 
of its data-points above the fitted lines are considered to be sin- 
gle young populations. 

4. Results 

First we stacked the candidates in three different, light- 
weighted stacks for the fitting; the total sample of "LAE best" 
and GALEX-detected non-AGN LAEs, as well as the total 
sample divided into a "red" and a "blue" sample (with Vj— i + 
colours greater or less than zero respectively). These stacks are 
similar to those in Nilsson et al. (2009a), but updated with spec- 
troscopic results and new photometry from the UKIRT J and 
CFHT K s images. The magnitudes for the three samples can 
be found in Table [5] Secondly, we also fit all individual ob- 
jects with at least one detection in the Spitzer bands and/or in 
the K s band. This sample consists of a total of 58 objects (8 
with only Spitzer detections, 19 with only a K s detection and 
31 with detections in both Spitzer and K s ). Redshifts are set to 
2.25 for all galaxies since varying the redshift within the range 
of the narrow-band filter does not change the fitting results sig- 
nificantly. 

For the main results presented here, the Monte Carlo 
Markov Chain SED fitting code was run with 20 000 itera- 



20 




1.0 1.5 
Old Age (Gyrs) 



Fig. 2. Mass fractions of the young populations at which an old 
population is discernible from a young population. Galaxies 
with mass fractions fit below these lines have two observable 
populations, for those with mass fractions fit above the lines 
it is impossible to discern any older population. The limiting 
mass fraction where the older population becomes visible are 
shown with different colours, depending on the ages of the 
young population. At very young ages of the young popula- 
tion, all older populations are invisible. The solid vertical line 
marks the age limit of the young population in our fits. The 
solid horizontal lines are linear fits to the limiting mass frac- 
tion dependencies on the older population age. The peak at age 
~ 1.3 Gyrs is a numerical artifact due to the limited age grid in 
the GALAXEV models. This does not affect our results as we 
use the solid, fitted lines for the measurements. 



Table 5. Magnitudes for fitting of the stacked SEDs. 





Total 


"Red" 


"Blue" 


# 


122 


80 


42 


Bj 


25.29±0.05 


25.22±0.04 


25.36±0.12 


v.i 


25.02±0.06 


24.92±0.11 


25.22± 0.09 


i+ 


25.00±0.14 


24.79±0.13 


25.51±0.17 


z+ 


24.92±0.22 


24.61±0.15 


25.46±0.23 


i 


24.02±0.24 


23.63±0.22 


24.47±0.45 


K s 


23.52±0.10 


23.18±0.09 


24.46±0.25 


Chi 


> 23.46 


23.38±0.12 


> 23.46 


Ch2 


> 23.46 


23.22±0.09 


> 23.46 


Ch3 


> 20.70 


22.79±0.13 


> 20.70 


Ch4 


> 20.65 


22.75±0.19 


> 20.65 



First row gives number of galaxies in stack. Errors are the largest 
of either the background variance noise or the result of bootstrap- 
ping error calculation. 



tions, and fitted the data with two single stellar populations. 
To find the probabilities in the parameter space, the first 2000 
runs were omitted to allow the x 2 to level out. The best fit pa- 
rameters are calculated as the first and second moments of the 
distribution in each dimension, weighted by the \r- The full 
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set of results of this run are found in Table [6] When interpret- 
ing the dust parameter, Ay, it is important to keep in mind that 
this is a galaxy averaged dust absorption. Different parts of the 
galaxy will be observable at different wavelengths. The results 
for the dust parameter will be dominated by the absorption of 
the restframe UV part of the spectrum, see also section 14.3.1 1 

When doing SED fitting, many of the fitted parameters can 
be degenerate with other parameters. In Fig. [3] we show, as 
an example, the contours of probability for the parameters of 
old population age, mass fraction in the young population, dust 
Ay, stellar mass, and young population age for two example 
objects (LAE_COSMOS_37 and 154), one that has a single 
young population and one that has a detectable older popula- 
tion (see also section l4~3l ). Of the parameters fitted for, the stel- 
lar mass is the best constrained. Secondly, the dust Ay is also 
well fit. Age of the old population, mass fractions and metal- 
licities are the least constrained parameters. As can be seen in 
Fig. [3] our calculation of the best fit value and its error bars 
is a good estimate of the distribution of parameters. No clear 
dependencies are seen. We will hereafter not discuss metallic - 
ities further, as they are in almost all cases completely uncon- 
strained. 

4. 1. Stack results 

The results for the stacked samples are shown in Table|6](top). 
A general comment for all results, both stacked and individ- 
ual, is that the old population age is very hard to constrain 
and the error bars are typically several hundred Myrs, also in 
cases where there is a confirmed older population. There is 
a large consistency between the stack results. The old/young 
ages are similar between the stacks. As expected, the red ver- 
sus blue stacks display diverging properties in dust and mass, 
with the red stack being dustier and more massive (Ay = 
1.04, M* = 7.3 x 10 9 M ) than the blue stack (Ay = 0.11, 
M* = 3.9 x 10 8 M ). The total stack has values in between 
those of the red and blue stack results. The results of these three 
stacks are also well fitted, with \r < 4. 

4.2. Interpretation of stack results 

When fitting SEDs at higher redshift, typically many objects 
are stacked in order to increase the signal-to-noise of the pho- 
tometric measurements. But the question remains if the proper- 
ties found fitting such a stack will be the average properties of 
the individual objects. To test this with our data, the 40 best fit 
objects were stacked, as well as the two sub-samples of single 
young population galaxies and those with an older population 
(see the sec. 14.3b . and these stacks were fitted in an identical 
way as for the individual objects. The three stack were fitted 
with best \r = [3.78, 4.73, 4.41]. The results are presented in 
Table 

The values from the mean of the individual fits and those 
from the stack are not in great agreement as seen in Table|7] The 
stellar masses give robust results in all cases, but the ages are 
about a factor of two younger in the stacks than the mean of the 
individual fits for the total sample and the old population sub- 



Table 7. Comparison of stack and individual fitting results. 





Age Ay 


log M* 




(Gyrs) (mag) 


CM©) 


All individually fit 


Mean of indiv. 


1.02 ±0.11 0.76 ±0.04 


10.39 ± 0.04 


Stack 


n 44 +o.40 nl O+0.17 

u -^-0.15 U ' iz -0.08 


10 23+ 09 


CTdiff 


1.40 3.66 


1.69 


Old pop. observed 


Mean of indiv. 


1.10 ±0.13 0.63 ±0.04 


10.41 ±0.04 


Stack 


0.56 ±0.16 0.14±g;J| 


10 42+ 06 


CTdiff 


2.62 3.37 


0.13 


Young pop. only observed 


Mean of indiv. 


0.012 ±0.001 1.51 ±0.04 


10.30 ±0.05 


Stack 


021 ±+ 042 1 67± +0 13 


10-45±g : gS 


°diff 


2.18 1.37 


2.37 



Table gives the best fit parameters for the mean and standard devi- 
ation of the mean of individually fit objects and the x 2 -weighted 
parameters of the light stack of the same objects. The results at 
the top of the table are for the whole sample of 40 best fit objects, 
whereas the bottom gives the results for the two sub-samples of 
galaxies with only a single young population observed (6) or with 
an older population observed (34). The trends of ages and Ays 
decreasing when fitting a stack are obvious in the total and old 
sample. 



sample. Further, the dust extinction is lower by about a factor 
of five in the stacks than in the mean of the individual fits. We 
interpret these results as being due to a smearing of the UV 
slopes of the SEDs of the galaxies in the stacking. It is clear 
that results from stacking of objects, especially those of the 
dust extinction, should be interpreted with caution and may be 
considered lower limits. 

4.3. Individual results 

The the best fit parameter values for all individually fitted ob- 
jects are found in Table [6] (bottom). In this section we remark 
on the general properties of all these objects, and on one special 
case. Of the, in total 58, objects fitted, four are detected in the 
GALEX data, of which three are spectroscopically confirmed. 
In Fig.|4]the histogram of the best fit \r are shown. An obvious 
break is seen at \r ~ 10 in the distribution. In Fig.[5j three ex- 
amples of SEDs and best fit spectra with very different xf- are 
shown. As explained in sec. [3] the xf- i s not expected to be unity 
and it is clear that fits with Xr ~ 10 st iH appear to be very good 
fits, as seen in the example of LAE_COSMOS_183, Fig. El We 
thus choose to define objects as "well fitted" if its calculated 
merit value (xl) is less than ten. Of the 58 fitted candidates, 40 
have x 2 r < 10- Two of these candidates are GALEX-detected, 
both of which are spectroscopically confirmed. For the remain- 
ing candidates with very bad Xr> about half have Spitzer data 
points lying significantly off from the other SED points (see 
e.g. LAE_COSMOS_147, Fig. 0, and the other half have un- 
usual SEDs, possibly consistent with lower redshift. Because 
of the reasons mentioned above, the results for the galaxies 
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Table 6. Full results of SED fitting run. 



Stacked samples 


Stack 


Old Age 
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log M« 
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Gyrs 
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Z/Z Q 


M 
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Total non-K s 
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1 in + 0.85 
1 -- L ' J -0.63 
i qi+0.75 
i - ol -0.77 

J - OD -0.81 


0.079™ 
0.018±g : gg2 

o.o46i;;»« 


0.43 +u ,-'$ 

U.OO_ 22 

0.38«;?° 


U.Ol_0.l8 

i n4+ 09 

1-U4_ 35 

11+ 014 

u " LJ --0.07 


u -*-0.2 

n 4+°-° 

u -*-0.2 

o 2+ 00 

u - -0.18 


io.o7±g:±§ 

9.87+g-?? 
q 5q+o.24 

a - oy -0.26 


4 67+ 1 ' 54 

4 - D '-0.52 

5 3fi+°' 83 






Individual results 


1 

10 
12 
15 
18 


76 +llb 

u - ' u -0.59 
U - ZO -0.04 

n ao+l-45 
u -°^-0.30 


074+ u ' U17 

U.U(^_ 030 

007+ 001 

U.UU(_ ooi 

n 022+ 005 

UiU "-0 .005 
u - u ' '-0.024 

017+ 007 

u - ul ' -0.005 


U ' JO -0.38 

> 0.26 

> 0.52 
04+ 016 
24 +0 ' 62 

u - zq; -0.20 


46 +U B::l 

n o n +0.04 
z ' ou -0.04 

1 76+ 011 

L - '"-0.10 

12 +01S 

u ' ±z -0.09 
9 -I q+0.16 
z - ±J -0.19 


0-4±»;l 
0.0051°;°, 
0.005t°;« 

1 0+°° 
0.005l°° 16 


10.161S:m 

10.641°;™ 
10.601°;^ 
10.74l°;°« 


*- 04 -0.81 
9q qf?+0.98 
^o.ou_ 60 

q Q4+2-02 

c: 79+I.28 
d - ' z -0.66 

32.001JJ 7 


14.11 ± 0.89 

77.51 ±5.16 
11.64 ±0.59 


G 



Columns show Xr- we ight ec l best fit parameters. Ages are in Gyrs and the mass fraction is the fraction of the total mass in the young 
population. Metallicity is fitted in steps, see the text for a description. The Xr is the reduced chi-square. The eighth column gives the dust 
corrected UV star formation rates for the candidates with Xr < 10, see also sec. 14.41 The final columns marks the entry with a G if the 
object is detected in the GALEX images. Full table is published as online material. 
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1000 10000 1000 10000 1000 10000 

Kesl (^) Kesl (A) K esl (A) 



Fig. 5. Illustration of the goodness of fit with various \r from three candidates. The Xr are 1, 8 and 95 respectively from left to 
right in the figure. As can be seen, also fits with \r ~ 10 are reasonably good fits. 



with bad Xr are not reliable and we will concentrate on the 
40 candidates with Xr < 10 in the further analysis. 

For all galaxies, the test shown in Fig. [2] was performed 
in order to determine if a galaxy was best described by two 
stellar populations, or could be best fitted with a single young 
population of stars. Hereafter, the "Age" of a galaxy is consid- 
ered to be the age of the young population, if no older popula- 
tion can be seen, and the age of the older population, if this is 
considered observed (cf. Fig.0. The results from the SED fits 
can be found in Fig. [6] Of the galaxies ~ 15% exhibit single 
young populations (15l|%). Conversely, ~ 85% can be seen 
to host at least one population of older stars. This is similar to 
the results of Guaita et al. (2010). If the stellar mass in the two 
sub-samples is summed up, 12% of the total stellar mass is in 
the young galaxies, and 88% is in the old. Hence, the stellar 
mass per galaxy is on average slightly higher in those galax- 
ies with an older population. The older ages of these multiple 
population galaxies range between 0.5 — 2 Gyr with a mean of 
1.10 ± 0.13 Gyrs, but are relatively unconstrained. The ages of 



the single, young population galaxies are 12 ± 1 Myrs. The dust 
Ay values for the whole sample range between 0.0 — 2.5 mag, 
with mean Ay = 0.63 ± 0.04 mag and 1.51 ± 0.05 mag for the 
old and young populations. The stellar masses are in the range 
8.5 < logM* < 11 M , with, for the old/young populations, 
logM* = 10.4 ± 0.04 / 10.3 ± 0.05 M . In Tabled the lower 
limits are given on the mass fractions in the young population 
in those galaxies with a single young population observed. 

4.3.1. LAE ULIRGs 

It was shown in Nilsson & M0ller (2009) that some of these 
candidates have infrared fluxes consistent with being ultra lu- 
minous infrared galaxies (ULIRGs). We now ask the question 
if these galaxies are located in a certain part of the parameter 
space fitted, or if they are evenly distributed. Of the non-AGN 
ULIRGs in the sample, four were fitted with good Xr an d these 
are marked with a red ring in Fig. [6] They appear evenly dis- 
tributed throughout the diagrams. Further, they are generally 
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above the average in stellar mass and all contain a significant 
older stellar population, but are very evenly distributed in the 
extinction parameter. The extinction in the UV/optical (Ay) 
was derived from the ratio of UV to infrared fluxes in Nilsson 
& M0ller (2009), and ranged between 4.8 — 7.2 magnitudes 
for the four that have been fitted here. Clearly, the two meth- 
ods, of determining the dust extinction primarily based on the 
UV slope or including infrared data points, reveal vastly differ- 
ent extinction values. The most likely explanation for this dis- 
crepancy is that they are due to dust geometry effects, as both 
methods depend differently on the assumption that the dust is 
distributed in a homogeneous screen in front of each object. 

4.3.2. LAE_COSMOS_94 

One very interesting example of the LAE ULIRGs is 
LAE_COSMOS_94, fitted with a best xl = 5.8. This object, 
which is confirmed to be a Lya emitter by the spectroscopy, is 
the one in our survey with by far the highest equivalent width 
(EW = 768 A). It has a Lya magnitude of 24.33 ± 0.10 
but is very faint in the optical broad-bands with magnitudes 
~ 26 and is even undetected in the z + and J bands, but then 
shows a very red SED in the K s and Spitzer bands, see Fig. [7] 
The best fit properties of this galaxy are age = 1.32^0 73 Gyr, 
A y = 2.38 ± 0.13, and log M* = 10.91±g;§Sj M . These 
properties are very similar to the average properties of all the 
candidates, except it has the largest Ay of all. 

Very high Lya EWs have been proposed to be signs of 
Population III stars (Schaerer 2003, Tumlinson et al. 2003), 
whereas LAE_COSMOS_94 is obviously not a Population 
III object. Very high EWs could also appear in dusty ob- 
jects, where the dust geometry allows Lya photons to escape 
(cf. Finkelstein et al. 2008). This object also appears to be quite 
extended in the narrow-band image, with a PSF-subtracted 
FWHM of 2'.'04, corresponding to 16 kpc at z = 2.25. The 
PSF-subtracted FWHM in the Bj/r+ bands are l'.'28/l'/04, or 
10/8 kpc. The narrow -band, i.e. Lya emission, is clearly more 
extended than the broad-band component, as has previously 
been found for other LAEs at both low redshift (Ostlin et 
al. 2009) and at high redshift (M0ller & Warren 1998, Fynbo 
et al. 2001, Finkelstein et al. 2010). This is further evidence 
that the high EW in this object is driven by resonant scattering, 
as the Lya emission is expected to be more extended than the 
continuum emission in those cases (Ostlin et al. 2009). 

4.3.3. Correlations in the data 

To test if any correlations exist between the typical parameters 
stellar mass, dust Ay and age in the full sample, a Spearman 
rank test was performed. The Spearman rank test allows to test 
the existence of any monotonous relationships in the data. The 
results of the Spearman rank test can be found in Table [8] No 
significant correlations are seen. 

With the data at hand, it is also possible to study depen- 
dencies between the Lya fluxes and equivalent widths and the 
parameters in the SED fit. In Fig. [8j the Lya fluxes and EWs 
of the galaxies are shown as function of stellar masses and dust 



Table 8. Results for the Spearman rank test on several parame- 
ter correlations. 



Parameters 


Spearman rank p 


Correlation prob. 


mass — Ay 


0.17 


86% 


age - Ay 


0.03 


68% 


mass — age 


-0.10 


86% 


Lya Flux — mass 


-0.20 


86% 


Lya Flux — Ay 


-0.23 


91.4% 


Lya EW — mass 


-0.08 


68% 


Lya EW - A v 


-0.16 


78% 



The Spearman rank number is a measure of the correlation in the 
data. Points along a perfectly monotonic function would have a 
Spearman rank of ±1. A positive rank number means a posi- 
tive correlation between the parameters, and a negative number 
an anti-correlation between the parameters. 



Ay. In the top panels in this figure, the Lya flux is shown as a 
function of stellar mass and dust parameter. A lack of objects 
with large Lya flux and large dust Ay is seen in the top right 
panel, although it is not statistically significant. In the lower 
panels of Fig. [8] the Lya restframe EW is shown as a function 
of stellar mass and dust. The object with extremely high EW 
is LAE_COSMOS_94, see also sec. l4~12l No correlation be- 
tween stellar mass and EW is seen, as proposed in Pentericci 
et al. (2009). The Spearman rank test was applied also to these 
data-sets, revealing only a weakly significant correlation be- 
tween Lya flux and Ay, see Table[8] 

4.4. Dust corrected star formation rates 

When the dust absorption for each galaxy is determined, the 
UV spectrum can be corrected and the intrinsic UV flux re- 
trieved. From the intrinsic fluxes, dust corrected star formation 
rates can be derived according to the conversion of Kennicutt 
(1998). The uncorrected SFRs were presented in Nilsson et al. 
(2009a), where the UV-derived SFRs of all candidates were 
found to be in the range 1—40 M Q yr . For the galaxies with 
good fits to the SED, the Vj band (used for the SFR calcula- 
tion) can be corrected with the best fitted Ay. Within this sub- 
sample of 40 galaxies, the maximum SFR^y before correction 
was 15 M Q yr -1 and the median ~ 8 M Q yr" 1 . After correc- 
tion, the range of SFRs in this sub-sample is 4 — 120 M© yr" 1 , 
with a median of 16 M Q yr" 1 . These values are more similar 
to z « 2 galaxies (such as BX/BM, BzK and DRG galaxies), 
although still in the lower two-thirds of the range covered by 
those galaxies (Reddy et al. 2005). The dust corrected star for- 
mation rates of the well fit objects are given in Table [6] 

5. Redshift evolution in properties 

We now wish to study how the properties of LAEs have 
changed with redshift, but results from SED fits of stacked 
photometry should be interpreted with caution, as noted in 
sec. 14.21 Stellar masses appear to be robust to stacking, but in 
this case there is a bias between stacks and individual fits in that 
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Fig. 8. Dependencies between Lya flux (upper panels) and EW (lower panels) and best fitted stellar masses (left panels) and dust 
Ay (right panels). Typical error bars are marked in grey. ULIRGs are marked with a red ring. The very high EW object in the 
lower panels is LAE_COSMOS_94, see sec. 14.3.21 No clear dependencies are seen, besides an apparent lack of objects at large 
Lya fluxes and large extinctions. 



individual fits are typically only done on the brighter galax- 
ies, which will then be biased towards more massive objects. 
Hence, to compare stellar masses, it is more advantageous to 
study stacked results, as these will be less sensitive to flux lim- 
itations. Further, even though the age of a stack was less than 
the mean age from the stacked objects fitted individually, the 
results were still relatively consistent with those from fitting in- 
dividual objects. In Fig. [9] we show the stellar masses and ages 
for stacks reported in several previous papers. In this figure, a 
clear correlation between age and stellar mass is seen, as con- 
firmed by a Spearman rank test assigning a 99.94% probability 
to it. It also appears the LAEs typically get older and more mas- 
sive with lower redshift. This trend will not be affected by the 
potential bias in younger ages found when stacking, as such a 
bias would be likely identical at all redshifts, and would thus 
only shift the relation in the age direction. Note though that 
there is an intrinsic degeneracy between ages and masses due 
to increasing mass-to-light ratios in older galaxies. 

We further see that the z = 2.3 LAEs are among the most 
massive LAEs so far found. To rule out a possible bias result- 



ing from continuum detection levels, we search the literature, 
from where the data-points in Fig.[9]are taken, for information. 
For each redshift, the continuum limit as close to restframe 
~ 4100 A as possible was noted. This corresponds to a lo- 
cation just redward of the Balmer break which should ensure a 
less dust-sensitive/more mass-sensitive data-point. The stacked 
continuum limit was then calculated from the 3a limit in the 
images used in each case, corrected for the number of objects 
in each stack. Finally, these limits were corrected for redshift 
dimming by converting to absolute magnitudes. Comparing 
these absolute, stacked continuum limits for the eleven mea- 
surements in Fig. |9]reveals that the stacked limit for the sam- 
ple here is the faintest of all, with an absolute stacked mag- 
nitude limit of -18.1 (3cr, AB) and with the other ten results 
having varying continuum limits between -18.2 — 20.7. Since 
the stellar mass in principle is directly proportional to the flux 
in the red part of the SED (excluding small dust effects), it is 
clear that any existing, smaller stellar mass objects would have 
been detected in our stack, and that the conclusion that z = 2.3 
LAEs are more massive than those at higher redshifts is robust. 



K.K. Nilsson et al,; The nature of z ~ 2.3 Lya emitters 



11 



Any bias, if present, would lead to lower average stellar masses 
derived for the z = 2.3 LAEs than for higher redshift LAEs, 
and therefore our conclusion of larger stellar masses at lower 
redshift is strengthened. 

To further study the apparent redshift dependence on the 
stellar mass and the ages, we plot the data-points from Fig. [9] 
separately for stellar mass and age as a function of redshift in 
Fig. [TO] The Spearman rank probability of a correlation in the 
data is 99.92% for age and redshift and 99.45% for stellar mass 
and redshift. The solid line in the figure represents linear fits to 
the data, with the following parameters. 

log Age (Gyrs) = 3.43 - 1.24 x z (3) 

logM* (M Q ) = 12.41 — 0.84 x z (4) 

In the interest of understanding any possible trends in the 
dust extinction with redshift we plot SED fitting results for the 
stellar mass and Ay in Fig.Qj] In this plot, only fits done on in- 
dividually detected objects are shown as it was found in sec. 14. 21 
that the dust extinction cannot be well determined with stacked 
fits. Unfortunately, this leaves only a few points in the plot, and 
no real trends are visible. 

Several difficult choices have to be made when preparing to 
fit SEDs of high redshift galaxies, besides the choice to stack 
or not, all of which may influence the final result (see also 
Gawiser et al. 2009). Firstly, the theoretical model must be cho- 
sen. We have here chosen the GALAXEV model, as this is most 
common in publications to date. Secondly, different groups fit 
different star formation histories (e.g. SSP, constant SFR, ex- 
ponentially declining SFR) which will influence the results. As 
described in Sect. 13.11 nebular emission makes a strong contri- 
bution to the SED at very young ages, but this effect is typically 
not implemented in SED fitting codes. Some publications have 
included the narrow-band flux in the fitting, in an attempt to 
model the dust extinction (e.g. Finkelstein et al. 2008), whereas 
most groups tend to either ignore the data-point, or use it to 
subtract the Lya flux from the broad-bands. Finally, one effect 
which may be the more difficult to circumvent, is that of dis- 
parate restframe wavelength coverage at different redshifts. At 
higher redshifts, the fits rely increasingly on the Spitzer data 
points to constrain the restframe optical and near-infrared. The 
uncertainties there lie in the low spatial resolution, and that, at 
fainter fluxes, the images quickly over-crowd. 

6. Discussion 

In this paper we have presented the results of fitting the SEDs 
of a set of Lya emitters at z = 2.25. The results presented here 
represent the single largest sample of individually studied Lya 
emitters to date. This is for the simple reason that galaxies at 
z ~ 2 are much more easily observable in the continuum than 
at higher redshifts. In Nilsson et al. (2009a), we proposed that 
LAEs at z ~ 2 are dustier and more massive than at higher red- 
shift, and that they contain more AGN. These suggestions are 
largely reinforced with the results presented here. With the new 
Chandra images, the number of detected AGN in the sample 
has increased to more than double what was found in Nilsson 



et al. (2009a). The LAEs in this sample have a typical Ay of 
~ 0.6 and a typical stellar mass of a few times 10 10 M Q . 

It is remarkable that almost all of the sample with a sin- 
gle young population have Ay > 1.0, whereas nearly all of 
the older objects have Ay < 1.0. This is a potentially inter- 
esting result, albeit potentially biased due to selection effects, 
e.g. that young low-dust galaxies will be faint in the bands red- 
ward of the Balmer break. Interestingly, no clear correlations 
were found between the Lya properties of the galaxies, and 
their stellar masses or dust contents. The apparent, though with 
low significance, anti-correlation between Lya flux and Ay, 
with fainter Lya fluxes at larger Ay, is a classic result, indicat- 
ing that the extinction in general affects the Lya emission in a 
strong negative way. 

It may be argued that the results presented here are biased to 
the brightest end of the galaxy sample, in that only objects in- 
dividually detected in the redder bands are fitted. For example, 
summing the average stellar mass found for the total stack, and 
comparing it to the sum of the individually fit objects Csec. l4.3b . 
we find that those 30% of the candidates contain 88% of the to- 
tal stellar mass among the galaxies. These 30% are obviously 
among the brightest, although the secondary criteria of hav- 
ing small Xr i s a l so imposed. However, the same bias exists in 
all other similar campaigns to understand the properties of high 
redshift galaxies. For this reason we also chose to fit the stacked 
fluxes, and tested if the fitting results were consistent between 
stacks and individual fits. Of the three main parameters re- 
ported here, the stellar mass was the most robust and can be 
trusted in all situations. Ages were slightly lower when stacks 
were fitted, compared to individual results, but the parameter 
that is most incorrectly determined for a stack is the dust ex- 
tinction. Due to this uncertainty, and the few results reported so 
far for individually fitted objects, it is difficult to draw any con- 
clusions regarding dust content evolution with redshift. As for 
the stellar masses, the stacked average masses found at z ~ 3 
are of the order 0.5 — 1.0 x 10 9 M Q , whereas the total stack 
here was best fitted with an average mass of ~ 1 x 10 10 M , 
about a factor of ten larger. We have also shown here that our 
stacked flux limit in the red continuum bands is comparable to, 
or in many cases deeper than, the stacked limits in previous sur- 
veys, and hence that the increase in stellar mass with redshift is 
robust. 

Further, we see no evidence for two populations of LAEs, 
with a red and a blue sub-sample. Instead, the results in Fig. [6] 
follow a smooth transition from little to a large amount of dust, 
and from smaller to larger masses. If a flux cut was made in our 
sample at a higher flux, only the most massive objects would 
be individually detected. The stack of the remaining objects 
would be dominated by smaller mass objects and so if this 
stack is also fitted for, it would appear as though there was a 
dichotomy of LAE classes. The smooth transition of properties 
seen in the sample here rather indicates that Lya emitters can 
be found in any type of galaxy at z = 2.25. Thus, Lya emitters 
are not a homogeneous type of high redshift galaxy, certainly 
less so than i.e. Lyman-Break galaxies (Nilsson et al. 2010). 
Selecting by Lya emission will, and does, find a cross-section 
of Lya emitting objects at each redshift and even though this 
technique is complementary to other high redshift galaxy se- 
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lection techniques, the most valuable comparisons are probably 
with other galaxies at the same redshift. For example, Shapley 
et al. (2005) found that LBGs at z ~ 2 had stellar masses 
log M* = 10.8 M Q , Ay = 0.86 mag and ages around 1 Gyr, 
with 15% having very young ages (< 100 Myrs). These values 
are more consistent with the results found for LAEs here, than 
comparing LAE results at different redshifts. To conclude, Lya 
emitters trace systematically different galaxies at different red- 
shifts, rather than comprising a certain class of galaxies, and 
can thus be used as a tool to find interesting objects at high 
redshift. 

7. Conclusion 

Based on the work presented here, studying the photometric 
properties of z = 2.25 LAEs, and fitting the SEDs of the ob- 
jects with spectral templates, the following conclusions can be 
made. 

- LAEs at z = 2.25 are generally more massive and possibly 
more dusty than at higher redshift, and their properties are 
quite diverse 

Using a sophisticated MCMC fitting method, the stellar proper- 
ties were derived for 40 well-fit objects, as well as three stacks. 
The LAEs were typically best fit with Ay = 0.0 — 2.5 mag 
and log M» = 8.5 — 11.0 M Q . We have shown that when mul- 
tiple populations of stars exist in the galaxies, the older pop- 
ulations are often outshone by the younger. Roughly 15% of 
the LAEs studied here are in such a strong star-bursting phase 
that any other older populations are invisible, whereas the re- 
maining 85% have "confirmed" older populations of stars. In 
section|5]we showed that LAEs at higher redshift (z > 3) tend 
to have smaller masses. 

- Larger AGN fraction among LAEs 

An analysis of the Chandra X-ray images available in 
COSMOS revealed 10 new AGN detections, increasing the 
AGN fraction in this LAE sample to 11%. 

- Stacking of objects does not reveal the average of the prop- 
erties of the individual objects 

SED fitting of high redshift galaxies is typically done on stacks 
of objects. We tested here the assumption that the stack reveals 
the average of the individual results. By stacking the individ- 
ually best fit 40 objects in different ways, we showed that the 
fitted stellar properties of this stack were not the average of the 
properties of the individual objects. In the process of stacking, 
the light from multiple populations of stars is smeared out, and 
particularly the Ay determination becomes highly uncertain. 

- Strong correlations between ages, stellar masses and red- 
shift were seen 

A correlation between the age and stellar masses with 99.94% 
probability was detected when comparing SED fitting results 
for LAEs at different redshifts. Further, a 99.92% probability 
correlation was found between age and redshift. This agrees 
well with a scenario where lower redshift LAEs are older galax- 
ies experiencing a second burst of star formation. 
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Fig. 3. Parameter dependencies between old/young popu- 
lation age, mass fraction, dust Ay and stellar mass for 
LAE.COSMOS.154 (left) and LAE.COSMOS_37 (right). 
Contours mark 68%, 95% and 99.7% probability. The green 
point marks the x 2 -weighted best fit parameter, with the la er- 
ror bar as found in Table [6] For the object in the left panels, 
the stellar mass is well constrained. The mass fraction in the 
young population is low, and hence the older population age is 
also well constrained. In this case, the younger population has 
a larger span of possible solutions, partially dependent on the 
Ay. This is an example of a galaxy for which the older popu- 
lation age is observed, as opposed to the object on the right. 
Here, the stellar mass and Ay are well constrained, but the 
young population age is very young, and the mass fraction in 
the young population is allowed to be any value above ~ 50%. 
Hence, the older population age is completely unconstrained. 
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Fig. 4. Histogram of Xr of non-AGN selected candidates. Right 
panel shows complete distribution with a bin size of 10, left 
panel a zoom on the range — 40 Xr with a bin size of 2. Filled 
histogram is for candidates that are not detected in the GALEX 
data. 



2.5 
2.0 

"1? 1.5 
o 
^E 

<? 1.0 
0.5 



o 
E 

o 

"53 

(7) 

en 
o 



o.o L_ 

0.0 



1 1.5 
1 1.0 
10.5 
10.0 

9.5 b 

9.0 

8.5 



r 

• 1 • 

• • •. • 
•* • . 

L . T~ 
r— * — i 



0.5 1.0 

Age (Gyr) 



1.5 



0.0 0.5 1.0 

Age (Gyr) 



1.5 



2.0 





® 

* • 

r 

• • 


• . ® 

• 


• 

• _ 
• 


- . •r ^ 


• 1 

• 




— • • 


• 



2.0 



1 1.5 










1 1.0 




• 




® - 


10.5 




• 

• 


• 

• 

• • 
• 

• • 


• — 


10.0 


• 
• 


; ••* ■ — 

• 


+* . 




9.5 




• 

• 


• 




9.0 




• 

• 






8.5 











0.0 0.5 1.0 1.5 2.0 2.5 

Ay (mag) 

Fig. 6. Results of 40 objects with xl < 10 (black points). The 
coloured data points represent the stacked SEDs, with green as 
the total sample, red and blue points are the red and blue sub- 
samples. Typical error bars are marked in grey. In the two top 
panels the error bars in the "Age" direction refer to those with 
ages above 100 Myrs. For the single young population objects, 
the typical error in age is ~ 15 Myrs. LAE ULIRGs (from 
Nilsson & M0ller 2009) are marked with a red ring. 
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Fig. 7. SED of LAE_COSMOS_94. It is very faint in the broad- 
bands, although with a large Lya flux, and a very red SED in 
the Spitzer bands. Red points are observed data-points, arrows 
mark upper limits. The solid line shows the best fit spectrum. 
The green data point is the narrow-band flux (Lya) for refer- 
ence. 
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Fig. 9. SED properties derived for stacked LAEs in several 
publications at different redshifts; at z ~ 3 with light blue 
points from Gawiser et al. (2006, 2007), Nilsson et al. (2007), 
Lai et al. (2008) and Ono et al. (2010a). Green points are 
at z = 3.7 — 4.5 from Finkelstein et al. (2007) and Ono et 
al. (2010a). The red point is at z ~ 5.7 from Ono et al. (2010b) 
and the brown point is from Ono et al. (2010b) for z = 6.5 
LAEs. The black star is the result of the fit of the total stack of 
LAEs (122 objects) at z = 2.25. Points without error bars were 
published without these values. 
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Fig. 10. Panels show the SED derived ages and masses as a 
function of redshift from the different surveys mentioned in the 
caption for Fig. [9] The solid lines are best fit linear relations. 





1 1.0 




10.5 










(n 


10.0 


cn 




D 




5 






9.5 


o 




~a> 




t/> 


9.0 


Qi 




O 






8.5 




8.0 



0.0 



0.5 



1.0 
Ay (mag) 



1.5 



2.0 



Fig. 11. Mean SED properties derived for individual LAEs in 
several publications at different redshifts; at z ~ 3/3.7 with a 
light blue/green point from Ono et al. (2010a) and Finkelstein 
et al. (2009a). Red points are at z ~ 5.7 from Lai et al. (2007), 
Pirzkal et al. (2007) and at z ~ 4.9 from Yuma et al. (2010) and 
the dark blue point is the mean of a sample of z ~ 0.3 LAEs 
from Finkelstein et al. (2009b), with confirmed AGN excluded 
(Finkelstein et al. 2009c). The black circles marks the mean in 
our sample of 40 best fit candidates; the filled circle is for the 
sub-sample with two populations and the open circle is for the 
single (young) population galaxies. 



